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ABSTRACT 

In this Letter we point out that redshift surveys can break the degeneracy between the 
galaxy bias, the power spectrum normalization, a^,o and the growth factor, without the 
need for external information by using a simple and rather general parametrization for 
the growth rate, the well known 7 parametrization and measuring the power spectrum 
at least at two different redshifts. We find that in next-generation surveys like Euclid, 
erg o an d 7 can be measured to within 1% and 5%, respectively, while the bias b(z) 
can be measured to within 1 — 2% in each of 14 equal-width redshift bins spanning 
0.7<z<2. 



1 INTRODUCTION 



The issue of constraining the cosmological parameters by 
employing the 3D galaxy power spectrum as a summary 
statistic of the matter density perturbations has been widely 
explored in the last few years. Indeed, the power spectrum 
estimated from a redshift survey can give a wealth of infor- 
mation through its shape, amplitude and radial anisotropy 
induced by peculiar velocity-driven redshift space distortions 
(RSD). In principle the shape can be used to constrain back- 
ground quantities such as the dark energy equation of state 
and the geometry of the system that are related to the ex- 
pansion history of the universe. The amplitude and RSD of 
galaxy clustering can constrain the growth of cosmological 
perturbations and the biasing relation between galaxies and 
mass, i.e. the relation between the spatial distribution of 
galaxy and the underlying mass density field. 

Constraining the growth is essential in order to discrim- 
inate dark energy models with the same background prop- 
erties but different physical origins. This makes the power 
spectrum a powerful tool in discriminating standard dark en- 
ergy models from modified gravity theories. However, sharp 
discrimination is possible only if galaxy bias can be esti- 
mated precisely or, at the very least, marginalized over effi- 
ciently. 



was realized that information contained in the redshift dis- 
tortions could tighten constraints over cosmological param- 



eters (Amendola et al. 2005 Wang 2008 Linder & Cahn 



2007) and RSD began to be considered as a standard, ad- 



ditional probe (Guzzo et al. [2008JI , able to constrain both 
the growth rate and the dark energy responsible for the ac- 
celerated expansion of the universe. The growth factor itself 
began to be regarded as useful information to be exploited 



rather than to be marginalized over (Amendola et al. 
|Wang et al^ [2010] ) 



2005 



However, as noticed in several works (Wang 2008 Per 



cival & White 2009 Song & Percival 2009 1 , the degenera- 



cies between the growth factor G(z), the (linear) bias b{z) 
and the power spectrum normalization (7g,o = °~s( z — 0) do 
not allow us to constrain all three quantities simultaneously. 
One can only measure combinations such as f(z)os(z) — 
f(z) G(z)a s , , b(z ) a s (z) = b(z)G(z)a 8 , or p(z ) = f(z)/b(z 
(e.g. |Ross et al " 



(20061; Guzzo et al. (20081; Blake et al 



Even though the work of Kaiser ( 1987 \ made it clear 



( |2011[ ) ) , where f(z) is the growth rate, related to the growth 
factor via f(z) = dlnG(z)/d\na and a = (1 + z)^ 1 is the 
scale factor. Th e growth rate was originally parametrized by 
|Peebles| {l980"t as f(z) = f7 m (z) 6 in a matter-dominated 
cosmology. This expression was later generalized as f(z) — 
f2 m (z) 7 to accurately describe the growth rate in a variety 
of cosmological models ranging from Dark Energy models 
to non-standard gravity theories ( Lahav et al.j 1991||Wang| 



that one of the most promising ways to determine the fluc- 
tuation growth is to exploit the redshift distortion effect 
(see Hamilton (19971 for a review), in the first pioneering 



I fc Steinhardt) [19981 [Amendola fc Quercellini| |2004| | kinder] 
|2005j [Polarski fc Gannouji 20081. This parametrization ob- 



works (Seo & Eisenstein 2003 20071 dealing with Fisher 



matrix formalism as a tool to constrain cosmological pa- 
rameters with redshift surveys, the redshift distortions were 
in fact considered a sort of "noise" disturbing the measure 
of the baryon acoustic oscillations peaks. The attitude was 
to marginalize over the RSD and the growth. Later on, it 



viously implies G(z) = J fi^/(l + 2) dz. 

Several authors have adopted different strategies to 
break, bypass or ignore the degeneracy between G(z), b(z) 
and (78,0 depending on the parameters they were interested 
in constraining. 

• In general one fixes (or assumes external priors for) one 
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of the parameters and then estimate the others. Different au- 
thors have fixed the clustering amplitude <rs,o to the value 
estimated from the cosmic microwave background (CMB) 



data to constrain f(z) and b{z) (e.g. Di Porto et al 



the growth parameter 7 (e.g. Di Porto et al.||2012 



et al. 



et al. 



2011 1 or both the bias and the growth (e.g. Di Porto 



2012 1 or 



Belloso 



2012). Instead of fixing ag : o to a single value, some au- 



thors have assumed CMB priors for this parameter to con- 
strain f(z) (e.g. Guzzo et al.|2008 l or 7 and b(z) (Gaztanaga 



et al. 



2011) 



One can assume that General relativity holds and fix 



7 to the A cold dark matter (ACDM) value (e.g. Seo & 
Eisenstein| [20031 



Wang et al. 



Amcndola et al. 2005; Ross et al 



2010 1. For example Hawkins et al. (20031 and 



2006; 



Cole et al. (2005) fix 7 to estimate the bias while Percival 



et al. ( 2004 1 assume a further prior on Q m to determine ag,o- 



• Alternatively, one can consider some measurable com- 
bination of the above parameters that can also discriminate 
models efficiently. A popular choice is the "mass weighted" 
growth rate f(z)ag(z) which provides a good test for dark 
energy models (e.g .|Song fc Percival|20"09"l| White et al.|2009| 
Percival fc White||2009| |Blake et al.||2011[ |Carbone et at] 



2011 



Wang et al. 12010 1 



It is a fact that galaxy bias, growth and clustering ampli- 
tude cannot be independently estimated from the observed 



P(k) alone (Percival & White 20091 and yet effective con- 
straints can be placed provided that these quantities could 
be accurately parametrized under fairly general hypotheses. 
Although the idea is somewhat implicit in some of the works 
quoted above, it seems to us it has not been clearly pointed 
out nor discussed thoroughly anywhere. The scope of this 
work is to make it explicit with a simple proof, that we 
present in Sec. [2] In Sec. [3] we present the Fisher matrix 
method that we employ in order to obtain the aforemen- 
tioned constraints on parameters. Results are presented in 
Sec. [4] and finally in Sec. [5] we draw our conclusions. 



However, let us notice that, since we are not interested 
in constraining D(z) and H(z) in every redshift bin, but 
only the parameters they depend on, these functions and 
their combination C(z) are not considered further in this 
Letter. 

Let us assume that the growth rate can be modelled as 
Slm(z) 1 . The growth index 7 needs not to be constant. The 
simple parametrization 7(2) = 70 + 712/(1 + z) has been 
shown to reproduce the time dependence of 7 in a variety of 
non standard gravity models ( Fu et al. 2009 1. For the sake of 



simplicity we consider in this Section the case of constant 7, 
but our conclusion remains valid also with time-dependent 
growth index, as we argue below. 

Cosmological parameters are determined by fitting 
Eq. |TJ to the observed galaxy power spectrum. From the 
spectral shape in the linear regime, Po, one can determine 
those parameters that describe the expansion history H(z) 
(and then D(z)) (h, f2 m ,o, Ode, wq, Wi, etc.). From the 
amplitude of the power spectrum one determines the combi- 
nation A(z) = G(z)B(z). If this analysis can be performed 
at two (or more) redshifts and the growth rate is modelled 
as / = fl m (z) 1 , then the degeneracy can be broken. Let us 
prove it for the simple case in which the power spectrum has 
been measured at two redshifts z± and z 2 . In this case, one 
determines A and R at two different epochs, Ai, Ri, A2, R 2 
and can solve the linear system 



A 2 
Ri 

R 2 



Gi( 7 )Bi 
G 2 {-y)B 2 

n m,l g 8,0 

Bi 

n m,2 g 8,0 

B 2 



^,i^8,oGi(7) 
A 1 

n^,a 8l0 G 2 ( 7 ) 



(2) 
(3) 

(4) 
(5) 



By assumption fi m , Zi has been already determined because 
it depends only on the background parameters. Then the 
ratio of the last two equations yields 



2 LIFTING PARAMETER DEGENERACY 

Let us consider a measurement of the galaxy power spectrum 
in redshift space at different epochs, i.e. in different redshift 
bins z. In the linear regime, the shape of the power spectrum 



and its redshift distortions can be modeled as follows (Seo 



& Eisenstein 2003 2007) 



P b s {z; k, n) 



G {z)b (2)ct 8i o 



Hz) 1 " . 



Po 



D 2 f(z)H(z) 
D*(z)H F (z) 
+P s (z) 

C(z)G 2 (z)B 2 (z) [l + Riz)^] 2 P + P s (z) (1) 



fi is the di- 
is the linear 



where B(z) = b(z)a a , , R(z) = 
rection cosine and Po = P(k, 2 = 0) 
power spectrum at the present epoch. The factor C(z) = 
D F {z)H(z)/(D 2 (z)H F (z)), where D(z) is the angular di- 
ameter distance and H(z) is the Hubble parameter, takes 
into account the difference in comoving volume between the 
fiducial cosmology - the one we use to convert observed red- 
shifts into distances - (subscript F) and the true cosmology. 
Finally, we model the shot noise contribution as an addi- 
tional factor P s (z). 



R1A1 
R 2 A 2 



Gi(7) 
G 2 (7) 



(6) 



where the only unknown is 7 and therefore we can solve for 
it. Then, from the above equations, one estimates ag,o and 
b{z) as 



08, = 



_A ± Rt_ 



and 



fi 7 

Km) = -^7- 



(7) 



We notice that if the power spectrum is measured in 
three redshift bins, then one can also constrain the time 
dependence of the growth index 7(2) = 70 + 712/(1 + 2). 
This procedure assumes that the size of the bins and the 
number densities of galaxies within are large enough to pro- 
vide effective constraints to avoid that parameters' degener- 
acy, although broken in principle, can reappear due to poor 
statistics. As we show in the next section, this is, however, 
not of concern for large surveys as Euclid. The proof that 
the degeneracy can be broken can be also obtained using 
the Fisher matrix formalism (Fisher, 1935; Tegmark et al. 
|1997[ ). If one assumes / = Q, m (z) 1 , a constant growth index 
and the model power spectrum of Eq. |l]) , then the deriva- 
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tives of the spectrum with respect to b(z), 7 and o"s,o 



9 In P obB 
dlnb Zi 
9 In Fobs 

9 In P bs 



= 2 



2 9JnG 2 r2 m (zj) 7 lnfl m (^j)^ 2 
97 b Zi + n m («i)T/x 2 

2 

O"8,0 



(8) 

(9) 
(10) 



are not degenerate because of their different dependence 
from z and Dropping the 7 parametrization, i.e. treat- 
ing G(z) as an extra free parameter for every bin, would 
not allow us to lift the parameter degeneracy and a singular 
Fisher matrix would result. We note that here we assumed 
the bias to be scale-independent which should be true on the 
large, linear scales we are focusing on here. A possible scale 
dependence would modify the shape of the spectrum and 
should be accounted for to avoid systematic errors in the 
estimate of the background cosmology parameters. Further- 
more, a scale-dependent bias could make the Fisher matrix 
degenerate again and cause the parameter degeneracies to 
return. Clearly, in order to break the degeneracy in presence 
of a bias scale dependence, one should parametrize it as, e.g., 
a simple power law of k, employing a number of parameters 
smaller than the number of distinct redshift bins. This issue 
will be studied in a future work. 



3 FISHER MATRIX FORECASTS 

We now apply the strategy outlined above to a specific mea- 
surement and adopt the Fisher matrix approach to evaluate 
the expected errors on the measured parameters. For this 
purpose we choose a reference cosmological model and run 



CAMB ( Lewis et al. 2000 1 to obtain the linear power spec- 
trum in real space Po- 

Then, the model power spectrum in Eq. [l] depends on 
a number of cosmological parameters, listed in Tab. [I] In 
the upper part we list parameters that do not depend on 
redshift. Redshift-dependent parameters are listed in the 
bottom part of the table. In our application we consider 
the more general case of a time-dependent growth index 
7(2) = 70 + 7i2/(l + z). We adopt a slightly unconven- 
tional fiducial ACDM model with 70 = 0.545, and 71 = 
(to be compared with the concordance ACDM values 
70 = 0.556, 71 = -0.018 ( |Fu et al.[ |2009[ )) and allow for a 
time-dependent equation of state for the Dark Energy with 
Wo — —0.95 and w\ = 0. 

In this work we do not consider a possible smearing 
of the wiggles associated with baryonic acoustic oscilla- 
tion (BAO) ( |Seo fc Eisenstein| [2007| . Further, we ignore 
non-linear effects that modify the spectral shape and the 



RSD pattern on small scales (e.g. Blake et al. (2011 1) by 
limiting our analysis to small wavenumbers k < fc max = 
min[fc cu t, fciin(z)], where fc cut is set to 0.2 h Mpc -1 . The 
value for fcim(-z) is set by requiring that the variance in 
cells a 2 (k\i n ,z) — 0.25 (for example, at z = 0.7 we have 
fci hl ~ 0.16/iMpc' 1 ). 

Having assumed our model power spectrum and set the 
parameter values of the fiducial model, we compute the ele- 
ments of the Fisher matrix as in|Tegmark ( 1997J) by integrat- 
ing over all modes below fc max . The result depends on the 



Table 1. List of the parameters that characterize the background 
cosmology and the model power spectrum. They are used in the 
Fisher matrix analysis. Redshift dependent parameters are listed 
in the lower part of the table. 



Redshift-independent 




fiducial values 


parameters 






Reduced total matter density 




0.271 ■ (0.703) 2 


Reduced baryon density 




0.045 ■ (0.703) 2 


Curvature density 


n k 





Hubble constant at present 


h 


0.703 


Primordial fluctuation slope 


n s 


0.966 


Dark energy eq. of state 




-0.95, 


Power spectrum normalization 


O"8,0 


0.809 


^-parameterization parameters 


70, 71 


0.545, 


Redshift-dependent parameters 






(in 14 z-bins) 






Shot noise 


Pa 









Derived from 


Bias 


log 6 


lOrsi et al.l 






polo} 



characteristics of the galaxy redshift catalogue used to com- 
pute the power spectrum. More precisely, it depends on the 
volume of the survey and on the galaxy redshift distribution 
dN/dz. In this work we take, as a reference case, the Euclid 



redshift survey as specified in the Red Book (Editorial Team 



et al. 2011 1 and at the website http://www.euclid-ec.org/. 
This survey will span a broad redshift range 0.7 < z < 2 that 
we split in 14 bins of Az — 0.1. The expected sky coverage 
is 15000 deg 2 . The expected dN/dz is given in the afore- 
mentioned Red Book, while the fiducial values for the bias in 
each redshift bin, for Euclid galaxies, are derived from |Orsi] 
etaLlpoTO). 



4 RESULTS 

Our analysis allows us to set simultaneous constraints on 
b(z), as and 7, when the power spectrum is measured at two 
or more redshifts. Tabs. [2p] display the expected la uncer- 
tainties on these (and other) parameters, listed in the first 
row. Errors on each parameter are obtained after marginal- 
izing over all other parameters in the table. Therefore they 
do not coincide with the length of the ellipses' axis shown 
in Fig. [l] which represent joint, rather than marginalized, 
probabilities. Table [2] shows the uncertainties on the red- 
shift independent parameters listed in the top row. In the 
second row the values refer to la uncertainties computed 
marginalizing over all other parameters while in the bottom 
row they are obtained through marginalization after fixing 
71 in order to compare our results to other, similar analysis. 
In Table[3]we list the expected uncertainties on the bias pa- 
rameter estimated at the different redshifts specified in the 
top row. The meaning of the second and third rows is the 
same as in Table [2] 

In the case of a time-dependent growth index, uncer- 
tainties on the parameters are remarkably small Aas = 
0.03 , A70 = 0.19 and b(z) = 2-4% depending on the 
redshift. They decrease even further if 7 is assumed to be 
constant (Aa 8 = 0.007, A7 = 0.03 and Ab(z) = 1 - 1.7%). 
In Fig. |T]) we plot the confidence regions of 70 and as- Con- 
tours refer to 68% and 95% joint probability levels. The blue, 
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dotted ellipses refer to the case in which we marginalize over 
all parameters, including b(z). The contour levels shrink con- 
siderably when one fixes the value of 71 (red, dashed) and 
become tiny in the case in which all cosmological parameters 
and the bias are fixed to their fiducial values (we leave only 
the shot noise free to vary). This last case represents the 
ideal situation in which the available complementary probes 
(e.g. SN la, CMB fluctuations, lensing etc) allow us to esti- 
mate all other parameters with very high precision. 

In Fig. [2] we show the analogous confidence contour lev- 
els in the as — b(z) (left-hand panel) and 7 — 6(2:) (right-hand 
panel) planes. As in the left-hand panel of Fig. [I] ellipses of 
different sizes refer to marginalization over different sets of 
parameters. The four sets of ellipses refer to different red- 
shift slices, as specified by the labels. We do not show all 
14 redshifts to avoid overcrowding. Keeping 7 constant and 
assuming a flat universe allows to reduce errors on as sig- 
nificantly but has little effect on galaxy bias uncertainties. 

Finally, we have tested the sensitivity of our result to 
fccut, i.e. to the cut we have imposed to exclude nonlinear ef- 
fects. Increasing fc cu t increases the number of fc-modes used 
in the analysis and reduces statistical errors. However, if k cut 
is too large, non-linear effects kick in and induce systematic 
errors. Decreasing fc cu t allows to apply linear theory but in- 
creases random errors. To find the best tradeoff one needs 
to estimate the power spectrum in some simulated galaxy 
catalog obtained from fully nonlinear N-body simulations. 
However, we can obtain an order of magnitude estimate for 
the relative importance of non-linear effects by re-computing 
errors at different values of fc cu t and check how much they 
change. 

The right-hand panel of Fig. [T] shows the result of 
such exercise in the 70 — erg plane. Ellipses filled with dif- 
ferent shades of blue represent the reference case of fc cu t = 
0.2 ft Mpc -1 (and correspond to the blue, shorth-dashed el- 
lipses in the left panel) . Pushing fc cu t up to 0.5 ft Mpc - (red, 
long-dashed) has a little effect, meaning that the decrease in 
statistical errors does not justify the risk of introducing sys- 
tematic errors driven by nonlinear effects. Reducing fc cu t to 
0.1 h Mpc -1 has a quite dramatic effect (purple, dotted el- 
lipses) . The large increase in the errors reflects the fact that 
with this fc cu t one cuts off all but one of the BAO wiggles. 
Finally, a value of fc cut = 0.15 h Mpc -1 (black, dot-dashed el- 
lipses) seems to represent a safer and yet acceptable option. 
However, the optimal choice of fc cu t depends on the under- 
lying cosmological model and on the accurate modeling of 
nonlinear effects. 



5 CONCLUSIONS 

The goal of this Letter is to point out explicitly that a two- 
point statistics like the power spectrum can provide inde- 
pendent constraints on the cosmological parameters o"8,o, 7 
and the galaxy bias, under the quite general assumption that 
the growth rate of density fluctuation can be parametrized 
as / = f2 m (z) 7 and the measurement is performed at two 
(or more) different redshifts. This result is not surprising. 
In fact it may appear even obvious, but to the best of our 
knowledge it has never been explicitly pointed out in the lit- 
erature. To assess how precisely these parameters can be de- 
termined after breaking the degeneracy we have performed a 



Fisher-Matrix analysis and explored the case of the next gen- 



eration Euclid redshift survey (Editorial Team et al. 2011 1 



We found that these quantities could be measured to the 
percent level. Constraints tighten considerably if a flat uni- 
verse is assumed or if 7 is taken to be time independent. 
We find that as and 7 can be measured to within 1%, 5%, 
respectively, while the bias b(z) can be measured to within 
1 — 2% (assuming flatness) in every redshift bin. 

This procedure can be applied to redshift surveys like 
WiggleZ ( Blake et al.||20l"l| that are deep enough to provide 
independent redshift shells in which to measure the power 
spectrum. Analyses of this dataset have provided estimates 
of combinations fas,oG(z) and /3(z) ( |Blake et al.| |2011[ ). 
We will investigate the possibility of using the WiggleZ to 
estimate as, 7 and b(z) in a future paper. 
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Par. 


h 




&b,0 


n 


^DE 


w 




70 


7i 


°"8,0 


la 


0.008 


0.004 


0.0008 


0.03 


0.009 


0.07 


0.31 


0.19 


0.49 


0.03 


1(7 


0.008 


0.004 


0.0008 


0.02 


0.009 


0.07 


0.28 


0.03 




0.007 



Table 2. Uncertainties on cosmological parameters, collectively indicated as Par. In the second row la values refer to 68% 
probability regions computed marginalizing over all other parameters while in the third row values are obtained through 
marginalization after fixing 71. 



2 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


la 


0.023 


0.022 


0.022 


0.023 


0.025 


0.027 


0.029 


0.031 


0.034 


0.038 


0.039 


0.041 


0.044 


0.046 


la 


0.010 


0.011 


0.011 


0.012 


0.013 


0.014 


0.016 


0.016 


0.018 


0.02 


0.021 


0.022 


0.024 


0.026 



Table 3. In this Table we show the la uncertainties over bias parameters in each redshift bin: in the second row la values refer to 
uncertainties computed marginalizing over all other parameters while in the third row values are obtained through marginalization 
after fixing 71 and Qk- 
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Figure 1. Left-hand panel: contours at 68% and 95% of probability for the parameters 70 and ergo when i) we marginalize 
over all other parameters (blue, dotted ellipses); ii) we marginalize over all other parameters after fixing 71 (red, long-dashed 
ellipses) and iii) we fix all other parameters. The last case is represented by the very small purple, solid ellipses nearly 
overlapping on the black dot which marks the fiducial model. The marginalized la errors for agfi are 0.03, 0.007, 0.002 for 
the 3 cases, respectively while those for 7 are 0.19, 0.03, 0.009. Right-hand panel: contours at 68% and 95% of probability 
for the parameters 70 and crg.o varying k cu t- 
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Figure 2. Left-hand panel: the three sets of ellipses show the contours at 68% and 95% of probability for the parameters and ergo 
and b(z) in 4 redshift bins when i) we marginalize over all other parameters (blue, dotted ellipses); ii) we marginalize over all other 
parameters after fixing 71 and !lg (red, long-dashed ellipse) and iii) we fix all other parameters (purple, solid ellipses). The black 
dots mark the fiducial values for the parameters and are centered at the central values of each redshift bin. Right-hand panel: as 
in left-hand panel, but for the parameters 70 and b(z). 



